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^ ' Abstract: Inverse nodal problem on Dirac operator is finding the parameters in the boundary 

D ' conditions, the number m and the potential function V in the Dirac equations by using a set of nodal 



points of a component of two component vector eigenfunctions as the given spectral data. In this study, 
we solve a stability problem using nodal set of vector eigenfunctions and show that space of all V 
functions is homeomorphic to the partition set of all space of asymptotically equivalent nodal sequences 
induced by an equivalence relation. Furthermore, we give a reconstruction formula for the potential 
function as a limit of a sequence of functions and associated nodal data of one component of vector 
eigenfunctions. Our method depends on the explicit asymptotic expressions of the nodal points and 
nodal lengths and basically, this method is similar to [T] and [5] which is given for Sturm-Liouville and 
Hill operators, respectively. 

^ : MSG 2000 : 34A55, 34L05, 34L20. 
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. 1. Introduction 

■ Inverse spectral problems have been an important research issue in mathematical physics. Using 

spectral data, different methods have been proposed for recovering coefficient functions in differen- 
tial equations [3], [3], [S], [B], [7j, [5], [S], [TD]. Normally, the spectral data contains the eigenvalues 
and a second set of data. In 1988, McLaughlin show that knowledge of the nodal points can de- 
termine the potential function of the Sturm-Liouville problem up to a constant [llj. This is so 
called inverse nodal problem. Numerical schemes were then given by Hald and McLaughlin [12] for 
reconstruction of the density function of a vibrating string, the elastic modulus of a vibrating rod, 
the potential function in the Sturm-Liouville problem and the impedance in impedance equation. 
Independently, Shen and his coworkers studied the relation between the nodal points and density 
function of string equation [T3]. Inverse nodal problems have been studied by many authors [14] , 

The Dirac equation is a modern presentation of the relativistic quantum mechanics of electrons 
intended to make new mathematical results accesible to a wider audience. It treats in some depth 
relativistic invariance of a quantum theory, self-adjointness and spectral theory, qualitative features 
of relativistic bound and scattering states and the external field problem in quantum electrody- 
namics, without neglecting the interpretational difficulties and limitations of the theory |16j . 

Inverse problems for Dirac system had been investigated by Moses [T7j, Prats and Toll |18j . 
Verde [19], Gasymov and Levitan [20], and Panakhov [21]. It is well known that two spectra 
uniquely determine the matrix valued potential function In [23], eigenfunction expansions 

for one dimensional Dirac operators describing the motion of a particle in quantum mechanics are 



investigated. In addition, inverse spectral problems for weighted Dirac systems are considered in 



For the Dirac system with a spectral parameter in the boundary conditions, one studied the 
properties of the eigenvalues and vector-valued eigenfunctions for the Dirac system with the same 
spectral parameter in the equations and the boundary conditions [25]. The sampling theory of 
signal analysis associated with Dirac systems when the eigenvalue parameter appears linearly in 
the boundary conditions was investigated [26] . One investigated a problem for the Dirac differential 
operators in the case where an eigenparameter not only appears in the differential equation but is 
also linearly contained in a boundary condition, and proved uniqueness theorems for the inverse 
spectral problem with known collection of eigenvalues and normalizing constants or two spectra 



Inverse nodal problems for Dirac systems had not been studied until the works of Yang, Huang 
[28j . They give reconstruction formulas for one dimensional Dirac operator by nodal datas. Later 
years, some problems were solved for inverse nodal problem of Dirac systems |16] . [29| . 

The organization of this paper is as follows: in section 1, we give some properties of Dirac 
system and quote some important results to use in main theorem. In section 2, we obtain some 
reconstruction formulas for potential function. Finally, we define do, dY,^^.^, Lipschitz stability and 
we show dsoj^ is a pseudometric on T,ijir. Then, we prove Theorem 3.1 in section 3. 

Consider the following Dirac system 

By'{x) + Q{x)y{x) = Xy{x),0 < x < vr, (1.1) 

with the boundary conditions 

(Acosa + ao) yi(0) + (Asina + 6o) 7/2(0) = 0, 

(Acos/3 + ai)yi(7r) + (Asin/3 + 6i)2/2(^) = 0, (1.2) 

where 

and y is a real valued, continuous function from the class Li[0, vr]. m,ak,bk{k = 0, l),a and /3 are 
real constants: moreover — § < a < ^ and — f < /3 < f [29] . 

The properties of the eigenvalues and eigenfunctions of the boundary value problem (1.1) were 
studied in [25]. Throughout their paper, Yang and Pivovarchik assume that 

oo sin a — bo cos a > 0, 

ai sin /3 - 61 cos /3 < 0. (1.4) 

Under this condition, the eigenvalues of the boundary value problem (1.1)-(1.2) are real and alge- 
braically simple |29] . 






Considering (1.3) in (1.1), we get 

V{x) + m 

V{x) — m 

and thus, equation (1.1) is equivalent to a system of two simultaneous first order differential equa- 
tions 

y2{x,X) + [V{x) + m]yi{x,X) - Xyi{x,X) = 0, 
-y[ix,X) + [Vix)-m]y2ix,X)-Xy2ix,X) = 0. (1.5) 

In generally, potential function of Dirac equation (1.1) has the following form 

Qix)-- 



Pn{x) Pi2{x) 

P2lix) P22{x) 



where Pik{x) (i, /c = 1, 2) are real functions which are defined and continuous on the interval [0, vr]. 
For the case in which ^12(2;) = P2i{x) = and pii{x) = V{x) + m, P22{x) = V{x) — m where m is 
the mass of the particle, the system (1.5) is known in relativistic quantum theory as a stationary 
one dimensional Dirac system or first canonical form of Dirac system [3] . 

We will use only one component yi{x, Xn) of the two component vector eigenfunctions y{x, A„) = 
[yi{x, Xn),y2{x, Xn)]^ . y{x,Xn) be the eigenfunction corresponding to the eigenvalue A„ of Dirac 
system (1.1)-(1.2). Suppose xtl^ are the nodal points of the component yi{x, A„) of the eigenfunction 
y{x,Xn)- In other words, yi{x'k^ , Xn) = 0. It is pointed out that the nodal points are repeated in 
accordance with the multiplicity of zeros. Let I^'' = Xn^, xi^^'^j , and denote the nodal length by 



We also define the function jn{x) to be the largest index j such that < xt^^ < x. Thus, j = jn{x) 



if and only if x G Xn , Xn ' for n > 0. 

Denote X = {xil^}. X is called the set of nodal points for the component yi(x,A„) of the 
eigenfunction y{x, A„) of the Dirac system. 

Now, we will give some results which is given in to use in main results. 

Lemma 1.1. [29] The spectrum of the problem (1.1)-(1.2) consists of eigenvalues {Anj^g^ 
which are all real and algebraically simple behave asymptotically as 

V c / 1 \ 

An = n-2 + - + - + — ,n^oo, (1.6) 
TT n yw^ J 

and 

V c ^ / 1 , 
A_„ = -nH \- U { ^ ] ,n ^ 00, 

TT n 

where 



TT 

m n\ ^0 sin a — bo cos a ai sin /3 — 61 cos (3 



TT 

V{x)dx + /3 - a,c= \ (sin 2a - sin 2/3) + 

2 27r 



TT TT 





(yi(x A) \ 
' , the solution of (1.1) satisfying the 



condition 

v(O.X) = I 

A cos a + oq 



!,(0,A)=(-;^='° " + '»'), (1.7) 



then, we have 

2/1 (x, A) = —A sin I Xx — V{t)dt + a 1 + — xcos I Xx — V{t)dt + a 

\ { / \ i J 

— mcosQsin Ax — / V{t)dt — oq sin Ax — / V{t)dt (1-8) 

, / \ / 



-bo cos Ax 



and 



y2{x, A) = —A cos I Xx — V{t)dt + a 1 H x sin I Xx — V{t)dt + a 

\ { J \ i , 

+m sin a sin I Ax - y V{t)dt | +aocosixx- J V{t)dt \ (1.9) 



/ \ 



-5o sin i^x-j y{t)dt J + O (^^^ , 



where r = \pX\ . 



Lemma 1.3. [29] For sufficiently large n > 0, the first component yi(x,A„) of the eigenfunc- 
tion y{x,Xn) of the Dirac system has exactly N{a,l3) nodes in the interval (0, vr) where 

n — 2, for a > and /? > or for a < and /3 < 
N{a, /3) = < n - 3, for a > and ^ < 
n — 1, for a < and /3 > 

Moreover, uniformly with respect to j G {1,2, N{a, (3)}, the nodal points and nodal lengths 
of the problem (1.1)-(1.2) satisfying the condition (1.7) has the following asymptotic formulas, 
respectively for sufficiently large n 



j7r+ y V{t)dt-a 



j i rn^xii^ m sin 2a ao sin a — 6o cos a r> ( ^ 

An 2A„ 2A„ A„ \X^ 



and 



i+i.i 



P'l = — + — 



m 



2jj,l 
I'll 



V{t)dt + —-^ + 



or by using asymptotic formula 



1 



+ 



v"^ — TT^C 



we conclude 



n-2 (n - 2fiT {n - 2)^7:'^ 
( 



V{t)dt - a 







V 



j-n + / V{t)dt -a\v 



( 



n-2 



(n-2)27r 



+ 



m^xi^ + m sin 2a + 2(ao sin q — 6o cos q) 
2(n-2)2 



+■ 



V 



m^i;^:^'"'^ + mf sin 2a + 2?j(ao sina — 6o cos a) ^ 1 



(n - 2)37r2 



(n - 2)37r2 



+ 



and 




[x-'n -xil ) 
2(n-2)2 



(n-2)37r2 ^'^U^ 
Now, we will consider (1.1) equation with following boundary conditions 

ni(0) cos a + U2(0) sin a = 0, 
ui(7r) cos /3 + ti2(vr) sin/3 = 0, 



(1.10) 



where < a,/3 < vr, and m is positive in (1.3). It is well known that the spectrum of the problem 
(1.1) with the boundary conditions (1.10) consists of the eigenvalues A„,n G Z which are all real 
and simple, and the sequence {A„,n G Z} satisfies the classical asymptotic form |28j 



An = n + - + — + 
vr n 



1 



(1.11) 



where 



>ci 



m(sin 2a — sin 2/3) + rr?T: 



f = /3 — a + y y(a;)(ix, i 



Let u{x,\) = (ui(x, A), ii2(x, A)) be the solution of the (1.1) with the initial conditions 

ni(0, A) = sin a, U2{0, A) = — cos a. (1-12) 
Then by the methods of successive approximations method [1] , there hold 



ui{x, A) = sin I Ax — y V{x)dx + a \ — y + O 



A2 



U2ix,X) = -cos\ Xx- I V{x)dx + a\ - — + [ — —\ (1.13) 



for large |A| , where [28] 

U (x, A) = —m sin I Xx — V{x)dx 1 cos a + ^-^x cos I Xx — V{x)dx + a 

\ i J ^ V { 

V{x, A) = msin I Xx — V{x)dx 1 sin a H xsin I Xx — V{x)dx + a 1 . (1-14) 

V { / V { / 

Lemma 1.4. [28j For suffuciently large \n\ , the i—th component Ui{x, Xn) of the eigenfunction 
u{x,Xn) of the problem (1.1), (1.12) has exactly |n| + l — i nodes in the interval (0,7r). Moreover, the 
asymptotic formulas for nodal points of first and second components of the eigenfunction u{x,Xn) 
as n ^ CO uniformly with respect to j £ Z are as following 

^" " A„ + A„ y ^^^^^^ A„ + 2A^ + 2A^ U^J ' ^ ^ ^ 

and 

^ . il^ ^fj n.).. - f + 1})!!^ , , o (JL) . ,1.1.) 

The asymptotic expansions of nodal lengths lil^ are the following for the problem (1.1), (1.12) 
by using (1-15) 



_^ 1 7 ... . , , {(-1)^+^ - (-l)nmsin2a , {{-ly+^xj^''' - {-ly xP^'W , ^ ^ ^ 
"A„ + A„ J ^^"^^^"^^ 2A2 + 2A,2 ^""Kxl 

In case of j = 2k{or j = 2k + 1), k €z Z; we obtain 

= i + i '7 '^'^'"^ ^ ^ ^ + o (3^) . (1.17) 



We can obtain /^'^ similarly as |?7-| — )■ oo by using definition of nodal lengths and (1.16). 
Theorem 1.1. Suppose that V G Li(0,7r). Then, for almost every x £ (0, vr), with j = jn{x), 



lim A„ I V{t)dt = V{x), 

n— >-oo J 



and 



where i = 1,2 and 



A. 



lim A„ / cos{2XnTTt)V {t)dt = 0, 

ra— >oo J 



n — 2, for the problem (1.1), (1-7), 
n, for the problem (1.1), (1-12). 



Proof: We prove this theorem by considering similar method in [T]. It is well known 
lim — I V{t)dt = V{x) and lim — / cos{2snt)V {t)dt = 

n— )'00 TT J n^oo vr J 

for almost every x £ (0, n) when n is sufficiently large |30j. If we use the fact 



we get 



lim — ^ = 1, 

n->oo TrAn 



lim — [ V{t)dt = lim A„ I V{t)dt. lim ^ 



lim A„ /" V{t)dt = V{x), 

n— >oo J 



and 



lim — / cos{2snt)V{t)dt = lim A„ / cos{2 Xn'JTt)V{t)dt lim — ^ 



lim A„ / cos(2A„7rt)y(t)dt = 0. 



The above equalities are still valid even if xil^ and oci^^'^ are replaced by X^^ and xl~^^'^. 

2. Reconstruction of Potential function by using nodal points 

In this section, we will derive some reconstruction formulas of potential functions for the prob- 
lems (1.1), (1.7) and (1.1), (1.12). Here and later, we denote (the grid lengths) to go with 
in the same way that (the nodal lengths) to go with [IJ. 



by 



m 



Theorem 2.1. Let V G Li[0,7r] be the potential function for Dirac system (1-1). Define Fn 

a) For the problem (1.1)-(1.7), 

fn-l 

Fn = {n-2)\Y, ri-2 

b) For the problem (1.1)-(1.12) 



k=l 



2{n-2) 



'n-l 



.A;=l 



n lb 



m 
2n 



^(n) ^ rr^^in) _ ^ 



Then, Fn converges to V pointwisely almost everywhere and also in Li sense. Furthermore, 
pointwise convergence holds for all the continuity points of V. 

Proof: 

a) We will consider the reconstruction formula for the potential functions for the problem 
(1.1)-(1.7). Observe that, by Lemma 1.3, we have 



1 



1 



and 



m 

2\r. 



vr 



V{x)dx + O 



Then, by using asymptotic expansions for eigenvalues below, we obtain 



m 

2\r, 



n-2 + 

n-2 + 
(n-2) 



/{'M n - 2 + O 
ih^ [n-2 



m 



m 



nl 2(n-2) + 0(i) 
2 /I' - 



vr 



2 n-2 \n 



vr 



li]^ {n-2- , 
" ^ 2 n-2 



vr 



+ o(l). 



Hence, to prove Theorem 2.1 (a), it suffices to show Theorem 2.2. 



(h) it can be proved similarly. To complete the proof of Theorem 2.1. (b), it suffices to show 
Theorem 2.3. 



An) 



Theorem 2.2. The potential function V G Li(0,7r) of the problem (1.1), (1.7) satisfies 

V{x) = hm lil^Xn - m^iTx 

for almost every x S (0, vr), with j = jn{x). 
Proof. By Lemma 1.3, 



V{x)dx 



2\i A 



vr xi, 
^ + — 



+ 0(t3 1, 



SO that 



,1 m'l^ 



2\i 



An - 71" 



An = An y y(x)dx + O (^-^ 



(2.1) 



We may assume Xn^ ^ x. By Theorem 1.1., if we take limit both sides of (2.1) as n — )• oo for 
almost X £ (0, vr), we get 

'•n 



V{x) = hm 



l]l An - m 



2 An 



vr 



An- 



Theorem 2.3. The potential function V G Li(0, vr) of the problem (1.1), (1.12) satisfies 

■^n ~r J'n j 



y(x) 



hm 

n— ^-oo 



2An 



vr 



An ± m sin 2a, 



for almost every x G (0, vr), with j = jn{x). 

Proof: It can be proved by similar process to Theorem 2.2. 
3. Main Results about Lipschitz Stability 

In this section, we solve a Lipschitz stability problem for Dirac operator. Lipschitz stability is 
about a continuity between two metric spaces. So, we have to first reconstruct to metric spaces. 
To show this continuity, we will use a homeomorphism between these spaces. Stability problems 
were studied by many authors [2], |31] . |32j . We will give a main theorem which execute that the 
inverse nodal problem for Dirac equation is stable with Lipschitz stability. 

We denote the space ^oir as a collection of all Dirac operators and the space as a collection 
of all admissible double sequences of nodes such that corresponding functions are convergent in Li. 
A pseudometric dso^r '^Dir will be defined. Essentialy, d^.j^^X-^-'-^) is so close to 

n-l 



do(X,X) = lim vr 



n 



2(n - 2) 



L 



k ) 



k=l 



where n > + 2 and L^"^ 
V2 



Yin) yin) _ -yin) -tfC") 



If we define X ^ X ii and only if dj^^^^^^X, X) = 0, then ~ is an equivalence relation on Scir 
and dsjjj^ would be a metric for the partition set ^^j^^ ~ ^Dir/ ■ Let SDin C Scir be the 
subspace of all asymptotically equivalent nodal sequences and let S^^^j,^ = ^Din/ ~ • Let $ be a 
homeomorphism the maps ^oir onto S^^^^. We will call $ a nodal map. 



(i) 



Definition 3.1. Let N' = N - {!}. We denote 

^Dir = {y £ ii(0,'7r) : V{x) is the potential function of the Dirac equation}, 
S Dir = The collection of the all double sequences defined as 

X = (X^ : A; = 1, 2, n; n G N', < X^ < X^ < ... < X^_^ < tt} , 
for each n G N. 

(ii) Let X GJ^Dir and define X = {XI'} where = X^^^ - X" and I^' = {XJ!,X]}^^) . We say 
X is quasinodal to some V G Li(0, tt) if X is an admissible sequence of nodes and satisfies 
(I) and (II) below: 

(I) X has the following asymptotics uniformly for /c, as n — >■ oo 



n-2 



+ 



n 



,k = l,2,...,n 



for the problem (1.1), (1.7). And the sequence of the functions 



'n-l 



.k=l 



n-2- 



m 



2(n-2) 



converges to y in Li. 

(II) For the problem (1.1), (1.12), X has the following asymptotics uniformly for fc, as n — >■ oo 



n 



X^ = — + Oi-],k = l,2,..., 



n 



And the sequence of the functions 



'n-l 



.fc=l 



converges to F in Li. 

Lemma 3.1. Let X,X eT,Di 
a) If X belongs to case I, then 



Til 



n zb 



m 
2^ 



TT 



n 



n-2 



+ 0( ^ ) ,fc = l,2,...,n 



If X belongs to case II, then 



Ll = - + o(-],k = l,2,...,n 
n \n ' 



h) Xn,k = \Xl-K\ = 0{l)- 

c) For all x G (0, vr), define Jn{x) = maks{k : Xj^ < x} so that k 
X G [Xj,Xj_^^]. Then, for sufficiently large n, 



Jn{x) if and only if 



\Jn{x) - Jn{x)\ < 1. 



Proof: 

a) For case I, we get 



rn Y» {k + l)TT kTT , ^ / 1 



n — 2 \n 



n — 2 n — 2 \n 
by using the definition of nodal lengths, and, similarly for the case (II), we obtain 

h) We will only consider case I. The other case is similar. By the asymptotic estimates 



\^k ~ ^k\ 



k 



klT kTT —r. 

+ --X, 



< 





kir 








+ 




n-2 




)+o 
















)■ 




n. 







n-2 n-2 



kiT 



n-2 



-XI 



= O 
= O 

c) Fix X G (0, tt). Let J = Jn{x) and J = Jn{x)- Since 

Jtt 



X:}<x< X]+, 



and 



X^<x<Xj^,^^ + o(-] =X^^Ux<X^^^ 



n-2 



n 



j+i " + n' ' 



when n is large enough, J + 1 > J and J +1 > J. Hence, — 1 < J — J < 1, then | J — j| < 1. 

Definition 3.2. Suppose that X,X E T^oir with and as their respective grid lengths. 



Let 



Sn{X,X)=TT 



n-2 



2{n-2) 



n-l 
k=l 



(3.1) 



Define 



do {X,X) = lim Sn {X,X) and d^^„. {X,X) = lim ^ • 

We obtain this metric by evaluating ||y — in Theorem 3.1. This definition was first made 
by [l]. Since the function f(x) = -^—^ — is monotonic, we have 

admitting that if do (X, X) = oo, then (is^.^^ (X, X) = 1. Conversely 

do {X,X)= d^oAX,^) 

We can easily proof this equality by using Law's method [I] , [2] . 

Lemma 3.2. Let X,X £ T^oir- 
O') ^Soir ^-^ ® pseudometric on Y^Dir- 

b) If X and X belong to different cases, then d^.^^^^. (X, X) = 1. 
Proof: It can be proved similar way with [1]. 

After following therem, we can say that the inverse nodal problem for the Dirac operator is 
Lipschitz stable. 

Theorem 3.1. The metric spaces (Jl/jiri II • II i) snd (^j^jj^.^^/ ~, (isj^j^) CLre homeomorphic to 
each other. Here ~ is the equivalence relation induced by dj:,^^^^. Furthermore 

,, _,, dsn iX,X) 



where d^^.^ {X,X) < 1. 

Proof: In view of Lemma 3.2. we only need to consider when X,X £ T,Dir belong to same 
case. Without loss of generality, let X, X belong to case I. In this case, we shall show that 

\\V-V\\^=do {X,X). 

According to the Theorem 2.1., Fn and Fn convergence to V and V, respectively. If we use the 
definition of norm in Li for the functions V and V, we have 



= {n-2) 



n-2 



2 

m 



2{n-2) 



j^in) _ jin) 

.Jn(x) Jn{x) 



dx + o{l). 



Hence by Fatou's Lemma, 

< (n-2) 

+ (n - 2) 



n-2 



m 



2(n-2)_ 



L 



n-2- 



m 



2(n-2)_ 



(n) _ jin) 
-H _ Yin) 

"M^) Mx) 



dx 



dx. 



Here, second and first term can be written as 



|y(n) y(n) 

\^Mx) ^J„(^:c) 



dx = o 



and 



I 



^Ju{x) ^Mx) 



k=l ^ ^ 



If we consider last equalities in (3.2), we get 



< (n-2) n-2 
+ (n-2) n-2 



m 



1 

2(n- 2)J ^ 

2 

m 



and 



n-2 



2(n-2) 



2(n- 2) 



^E|Li.-r^K„(i,) 

fc=i ^ ^ . 



n-l 



5;|L^-L^|+o(l). 



Similarly, we can get 



n- 2- 



m 



2(n-2) 

If we consider (3.3) and (3.4) together, we obtain 



n-l 



5^|L^-L^|+o(l). 



n-2 



fc=i 



n n-l 



(3.2) 



(3.3) 



(3.4) 



2(n-2)J^^ 

The proof is complete after taking the limit as n — >^ oo. 
3. Conclusion 

In this study, it was given a Lipschtiz stability of inverse nodal problem for Dirac operator 
by using zeros of one component of two component vector eigenfunctions. Especially, two metric 
spaces was defined and it was shown that these two metric spaces were homeomorphic to each 
other. These results are new and can be generalized. 
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